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ABSTRACT 

We solve the RXTE X-ray lightcurve of the extremely luminous and massive star 
i] Carinae with a colliding wind emission model to refine the ground-based orbital ele- 
ments. The sharp decline to X-ray minimum at the end of 1997 fixes the date of the 
last periastron passage at 1997.95 ± 0.05, not 1998.13 as derived from ground-based 
radial velocities. This helps resolve a discrepancy between the ground-based radial ve- 
locities and spatially-resolved velocity measures obtained by STIS. The X-ray data are 
consistent with a mass function f(M ) « 1.5, lower than the value f(M) « 7.5 previ- 
ously reported, so that the masses of ?/ Carinae and the companion are M n > 8OM 0 
and M c ~ 3OM 0 , respectively. In addition the X-ray data suggest that the mass loss 
rate from 17 Carinae is generally less than 3 X 10~ 4 M 0 yr _1 , about a factor of 5 lower 
than that derived from some observations in other wavebands. We could not match 
the duration of the X-ray minimum with any standard colliding wind model in which 
the wind is spherically symmetric and the mass loss rate is constant. However we show 
that we can match the variations around X-ray minimum if we include an increase of a 
factor of ~ 20 in the mass loss rate from r) Carinae for approximately 80 days following 
periastron. If real, this excess in M would be the first evidence of enhanced mass flow 
off the primary when the two stars are close (presumably driven by tidal interactions). 

Our interpretation of the X-ray data suggest that the A SC A and RXTE X-ray spectra 
neai the X-ray minimum are significantly contaminated by unresolved hard emission 
{E > 2 keV) from some other nearby source, probably associated with scattering of the 
colliding wind emission by circumstellar dust. Based on the X-ray fluxes the distance 
to 17 Carinae is 2300 pc with formal uncertainties of only ~ 10%. 
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1. Introduction 

Discovery of cyclical variations in the strength of the He I 10830 A line (Damineli 1996), the 
X-rav flux (Corcoran et al. 1995; Ishibashi et al. 1999) and radio emission (Duncan et al. 1995) 
spurred the suggestion that the supermassive star rj Carinae is in fact a binary system (Damineli 
1996; Damineli, Conti and Lopes 1997; Corcoran et al. 1997) in which the variable phenomena 
over all wavelengths is driven by wind-wind interactions. Since the formation and evolution of 
massive binaries proceeds via mechanisms distinct from those available to single stars (Vanbeveren, 
De Loore, and Van Rensbergen 1998), the presence of a companion would have a profound impact 
on our understanding of the formation and evolution of ?/ Carinae in particular and of extremely 
massive stars in general, and it is essential that we resolve this issue one way or another. 

However, the “binary model” as proposed faces three major observational difficulties. The 
first concerns the variation of the X-ray emission from the source as measured by the Rossi X-ray 
Timing Explorer ( RXTE , Bradt et al. 1993). First attempts to model the observed RXTE 2-10 keV 
X-ray variability (Ishibashi et al. 1999) showed qualitative similarities to the observed variations 
in flux and column density but could not describe the X-ray variability near the X-ray minimum 
if the orbital elements derived from ground-based radial velocity data (Damineli, Conti and Lopes 
1997; Davidson 1997) were assumed. A second problem is that the column density to the X-ray 
source as measured by the A SC A satellite (Tanaka, Inoue, and Holt 1994) is lower, not higher, 
than the value outside eclipse and no increase in flux was seen to energies of 10 keV, in contrast 
to model predictions (Corcoran et al. 2000). Finally, expected radial velocity variations (Damineli, 
Lopes and Conti 1999) were not confirmed by subsequent HST/STIS observations of the stellar 
emission resolved from contamination by nearby circumstellar emission (Davidson et al. 2000). A 
more recent solution of the ground-based data (Damineli et al. 2000) does better to reconcile the 
velocity curve with the STIS measurements but does not adequately resolve the problems with the 
observed X-ray variations. Since the only observable spectral features are emission lines, there is 
the real possibility that the variation of the emitting gas may not faithfully represent the motion 
of the star, which means that deriving a true stellar radial velocity curve from the emission line 
analysis alone may be difficult no matter what spatial resolution is available. 

A more recent analysis of the ATE data by Ishibashi (1999) showed that the match between 
tlie model and the X-ray lightcurve could be improved if the orbital elements were relaxed from the 
ground-based values. Here we try to use the X-ray lightcurve to attempt to refine the orbit. We 
find that we can generate an improved match to the RXTE fluxes if the time of periastron passage 
is earlier, the eccentricity higher and the mass function lower than given in the Damineli et al. 
(2000) orbital solution. Most importantly, we find that the width of the X-ray minimum requires 
a suppression of the observed flux for an extended period following periastron. One way to do this 
is by a temporary large increase (by a factor of 20) in the mass loss rate from i] Carinae for a brief 
period just after periastron passage. In the rest of the paper we describe our analysis of the X-ray 
lightcurve. In §2 we discuss the data extraction. In §3 we describe the colliding wind model and 
discuss our attempts to solve the X-ray lightcurve, and in §4 we present our conclusions. 
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2. Observations and Data Reduction 

We use data obtained by the RXTE Proportional Counter Array (PCA). Data extraction, 
and instrumental background correction, are as described in Ishibashi et al. (1999). Most of the 
early data were obtained using 3 proportional counter units (PCUs), though we also include scaled 
2-PCU mode data in the lightcurve, since most of the data available after 1998 is in 2-PCU mode. 
To convert these rates to flux, we use a linear conversion derived by comparison of RXTE and 
spatially resolved X-ray fluxes measured by contemporaneous observations of 7 / Carinae with the 
.1 SC A X-ray satellite. There are four ASCA observations which overlap RXTE pointings within 
2 days or less. We derived source fluxes from the ASCA data by direct fitting of the extracted 
spectrum corrected for instrumental background and for sky background (where “sky” is the region 
beyond the 3' radius circle used to extract the ASCA source spectrum). We then compared these 
fluxes to the RXTE count rate corrected for instrumental background. Figure 1 shows the ASCA 
fluxes and the RXTE net rates. A fit to these data yields flux = 0.26(A- 7.66), where flux is the 
2-10 keV flux seen by RXTE (in units of 10~ n ergs cm -2 s -1 ), and N is the RXTE 3-PCU rate 
corrected for instrumental background. Thus the sky background contributes about 7.66 counts 
s -1 to every 3-PCU observation. This value is slightly lower than the value used in Ishibashi et al. 
(1999) by about 1 count s -1 per 3 PCUs. 

We also included RXTE “quicklook” data for the period since mid 1998. For each observation, 
we extracted PCA rates for the intervals when the satellite was on target; we then corrected these 
data for instrumental background in the same manner as used for the processed data. Since many 
of these data were obtained in 2 PCU mode we performed an additional scaling to correct the 2 
PC'U data, to approximate 3 PCU rates. We compared these net rates to the RXTE fluxes derived 
from the processed data in the intervals of overlap. We found that the scaling from the corrected 
quicklook net rates to observed 2-10 keV flux is flux - 0.29(N gl - 8.86), where flux is the 2 10 
keV flux (in units of 10 -u ergs cm -2 s -1 ) and N q i the corrected quicklook rates. 

Figure 2 shows the RXTE lightcurve in the interval February 1996 through March 2000. This 
lightcurve includes the data published in Ishibashi et al. (1999) and more recent observations. For 
our purposes, the significant characteristics are these: 1) a gradual increase in count rate starting 
at 1997.0; 2) a rapid rise to maximum and even faster decline to minimum at the end of 1997; 
3) a gradual brightening of the X-ray flux coming out of the X-ray minimum; and 4) asymmetric 
maxima just prior to and just after the X-ray eclipse. Also apparent in the lightcurve are higher 
frequency “flares” (Corcoran et al. 1997; Ishibashi et al. 1999); in what follows we make no attempt 
to try to explain or fit these higher-frequency variations. 


3. X-ray Lightcurve Analysis 

The colliding wind emission model we use is based on the Usov (1992) model used in Ishibashi 
et al. (1999). We fix the period of the binary orbit at 5.52 years (Damineli 1996), since this value 
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scorns to describe the lie I 10830A variations (Damineli 1996), the radio variability (Duncan ct al. 
1995) and the recurrent X-ray low state (Ishibashi et al. 1999). We take the the orbital elements 
a-s given in Damineli et al. (2000) as a starting point, and, since the variations in He I 10S30A 
and the radio emission suggest an eclipse of a. source of UV radiation presumably provided by 
the companion (Damineli, Conti and Lopes 1997; Duncan et al. 1995), we consider solutions in 
which the companion is a strong source UV source, which restricts the mass of the companion to 
M c > 15M 0 . We further assume that the winds collide at terminal velocities, and take = 500 
km s _1 (Hillier et al. 2000) and V’oo^ = 2000 km s _1 as reasonable values for the terminal velocities 
of )j Carinae and the companion, respectively. 


3.1. A Constant M Model 

Table 1 lists the parameters we used to derive the model curves presented in figure 2. After 
some experimenting, we adopted mass loss rates M n = 1 x 10 -4 M@ yr _1 for ?? Carinae, and 
Me = 1 x 10~ 5 M @ yr -1 for the companion. We adopted 8OM 0 as the mass for ?? Carinae. We 
found that a companion mass of M c ~ 30 M 0 provided a suitable fit to both the X-ray lightcurve 
and the STIS radial velocities. These component masses are similar to those described by Davidson 
(1999). However, our analysis is not very sensitive to the values of the component masses, and we 
find that we can fit the X-ray and STIS data with larger values of q Carinae’s mass, or smaller 
component, masses, or some combination. We allowed the eccentricity e, time of periastron passage 
T 0 , and the terminal velocity of the wind of q Carinae I / 00i , y to vary to improve the fit to the RXTE 
lightcurve. We found that in order to match the rapid drop in X-ray brightness, the value of u> 
needed to be near 270° with periastron passage T 0 = 1997.95, with an estimated uncertainty in T 0 
of ±0.05 years. These parameters ensure that the line of sight at periastron passes through the 
thickest part of the wind from q Carinae if i < 60°. While our adopted value of u) is the same as that 
of Damineli et al. (2000), the sharp decline to minimum in the X-ray data fixes periastron passage 
to occur about 73 days earlier than the Damineli et al. (2000) value. We also derive an eccentricity 
c = 0.9, which is higher than the highest previously published value, e = 0.8 (Davidson 1999). The 
smooth solid curve in figure 2 is the colliding wind emission based on our adopted parameters; the 
model based on the Damineli et al. (2000) values given in table 1 is shown by the thick broken line. 
Our model is an improvement over those published previously (Pittard et al. 1998; Ishibashi et al. 
1999) in that the’ decline to minimum and the asymmetry of the shoulders of the eclipse are better 
matched (though the model badly fails to match the observed eclipse duration). The ratio of the 
model flux to the observed X-ray flux determines the distance to q Carinae as D ss 2300 pc. The 
formal uncertainty in the distance we derive is primarily influenced by uncertainties in fixing the 
X-ray flux level outside of the eclipse, and is only about 10%. However, uncertainties in the model 
parameters and other simplifications in the analytic description of the X-ray luminosity may make 
the real uncertainty in the distance somewhat larger. 

The values of M listed in table 1 are rather uncertain; in particular we have adopted a rather 
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large value for t he mass loss rate of the companion, i.e. M c = 1(T 5 M@ yr -1 , which would suggest 
t hat the companion is an evolved WR-type star. While the adopted values provide a fairly good 
match to most of the observed X-ray variability (outside of X-ray minimum), these mass-loss rates 
are not uniquely defined. Other combinations of M v and M c can also provide an adequate match 
to the data, especially if we allow some latitude in choice of terminal velocities. In order to see if 
any constiaints could be placed on mass loss rates, we ran models with differing values of Kf for 
bot h stars, with the wind terminal velocities and orbital elements held constant at the tabulated 
values. We found that the mass loss rate from r] Carinae is constrained to be M v < 3 x 10 - ' 1 
M® yr -1 in order to match the gradual rise to maximum in the 1997.0-1997.9 interval. Models 
with M v larger than this limit show significant absorption of the 2-10 keV X-ray emission months 
prior to periastron, which suppresses the rise in L x up to periastron. In addition these models 
generally yield columns of N H > 10 23 cm -2 at all orbital phases, while the observed colums are 
typically Nfj ~ 3 X 10 22 cm 2 for most of the orbit (?). As an example of this effect, we include in 
figure 2 a model with M n = 3 x 10" 4 M 0 yr" 1 . This model is basically flat through 1997 except 
for a very brief increase in observed flux very near periastron, in contrast to the observed X-ray 
variation during this interval. It is more difficult to put constraints on the mass loss rate from the 
companion, since the X-ray lightcurve in the 1996-2000 interval is dominated by the wind from tj 
Carinae. Observations near inferior conjunction/apastron (in 2002.25) may allow us to put tighter 
constraints on M c as the companion moves in front of 7/ Carinae. 


3.2. An Enhanced M Model 

While the new orbit model provides an adequate match to the decline to minimum, the width 
of the eclipse in the model is still much narrower than the observed X-ray eclipse. This is mainly due 
to the fact that, for any significant orbital eccentricity, the stars move very quickly near periastron 
so that the line of sight quickly passes through less dense regions of the wind from i] Carinae as the 
orientation to the colliding wind shock changes. As noted by Corcoran et al. (2000), however, the 
anomalously small column density measured by ASCA during the X-ray eclipse, and the lack of 
any enhancement in emission out to 10 keV at that time, may well indicate that the colliding wind 
source is completely hidden, and that the measured emission is simply contamination from another 
source of keV X-rays within the source extraction region (i.e. within 3' from 7; Carinae). Corcoran 
et al. (2000) noted that a column density Nj-i > 10 21 cm -2 would be necessary to completely hide 
the colliding wind source. A column so large might imply either a non-spherical wind or disk 
around 7/ Cannae (Ishibashi et al. 1999), or a mass loss increase near periastron if the wind from ?? 
Carinae is spherically symmetric away from the wind-wind collision shock. We attempted to model 
the X-ray lightcurve by allowing for a brief period of enhanced mass loss from 7 7 Carinae after 
periastron. Our best fit is shown as the solid smooth curve in figure 3. In this model we assume an 
inc lease of a factor of 20 in A7, ; in the interval T 0 < i < 77, -f 0.22 J years, i.e. we assume an interval 
of mass loss enhancement by a factor 20 starting at periastron passage and lasting for about 80 
days. The model shown in figure 3 provides a much better match to the observed duration of the 
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X-ray eclipse. For comparison we also calculate a similar “enhanced- A/,” model with the orbital 
elements of Damineli et al. (2000). This model is shown as the broken curve in figure 3. Figure 1 
shows the variation in the model N H relative to the constant M n model. 


3.3. Expected and Observed Radial Velocities 

From their STIS data, Davidson et al. (2000) determined velocity shifts in the range —23 < 
- < -4 km s -1 for the Pa 3-7, 3-8, 3-10, and 3-11 emission lines, where v 0 7 are the velocities 
they mea.su led at <0/6 of the peak of the emission line. The velocities uq .5 at 50% of the emission 
line peak lie in the range —3 < At>o.s < +4 km s _1 . The emission line velocity shifts are much less 
than the expected shift Ae = (u(1999.1) - v(1998.2)) = +35 km s" 1 based on the orbital solution 
of Damineli, Conti and Lopes (1997) and Damineli, Lopes and Conti (1999), which Davidson et 
a.1. (2000) used for comparison to their STIS data. The orbital solution in Damineli et al. (2000) 
implies that u(1998.2) = +26 km s _1 and t>(1999.1) = +24 km s" 1 or An = -2 km s~\ in better 
agreement with some of the STIS observations, though this solution still has serious problems in 
matching the observed X-ray lightcurve, as shown in figures 2 and 3. Figure 5 shows the Damineli et 
al. (2000) radial velocity curve, and the radial velocity curve based on the “X-ray orbital elements” 
in Table 1. The vertical lines show the range of v 0 , 7 from Davidson et al. (2000) for their Paschen 
line measurements. The v 0 5 velocities show an even larger scatter than the v Q _ T velocities. Both 
the Damineli et al. (2000) radial velocity curve and the curve based on the X-ray analysis show 
reasonable agreement with the STIS velocities. 


4. Conclusions 

Our attempts to model the RXTE lightcurve of i] Carinae strongly suggest that: 

1. t) Carinae is a. binary star in which the observed properties of the 2-10 keV X-ray emission 
are primarily due to wind-wind collisions; 

2. if the luminosity of the X-ray source varies as the inverse of the separation of the stars, then 
the duration of the X-ray minimum is very difficult to fit in any colliding wind model with 
significant eccentricity and i < 60° unless the observed flux is artificially suppressed for an 
interval of time following periastron, possibly due to an enhancement in mass loss rate at 
that time. 

3. the sharp decline to minimum requires that the periastron passage occurred at the end of 
1997, and not in early 1998. 

It could only be possible to keep the time of periastron passage as late as 1998.13 if there were 
a significant mass loss enhancement prior to periastron passage. Our preliminary exploration of 
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.such models suggests that they do not match the observed decline from X-ray maximum to X-ray 
mininum as well as the model presented here. It is of course harder to conjecture physical reasons 
foi an enhancement prior to periastron passage, though full illumination of this point requires 
detailed modeling of the tidal forces on 77 Carinae. The duration of the X-ray eclipse in the model 
could be extended without an enhancement in M n if the inclination is closer to 90°, since then the 
line of sight to the X-ray source passes through very dense regions of the wind from 77 Carinae for 
an extended period near conjunction. However large values of inclination seem ruled out by lack 
of optical variability. It is worth noting that, for values of i near 50°, and if the wind eclipse is 
complete, then conjunction must occur near periastron, since otherwise there is not enough wind 
material along the line of sight to produce much absorption at energies above 5 keV. 

In a very eccentric orbit, the increase in tidal force near periastron provides an obvious expla- 
nation for an increase in M v at that time, especially if the outer layers of 77 Carinae are not tightly 
bound to the stellar core. The importance of tidally driven mass loss from 77 Carinae has been 
considered previously (for example, Davidson 1997). Even if e = 0 . 9 , the periastron separation is 
still about 1.5 AU = 32O/? 0 , so the stars are not particularly close. However the radius of 77 Carinae 
is currently not well constrained, and probably lies in the range 70 - 450 R e (Davidson 1999), so 
even in the conservative case the separation is only about 5 times the stellar radius of 77 Carinae. 
It is interesting to speculate that some of the line profile variability observed near periastron is 
actually produced by a combination of a change in radial velocity and a change in circumstellar 
density; if so, then radial velocity curves may be biased in some way near periastron passage. 

Perhaps one of the more controversial results of our X-ray analysis is the rather low mass loss 
rate we derive for 77 Carinae. Typical values of M„ are generally ~ 10 _ 3 M@ yr -1 , not ~ 10 “ 4 A/ Q 
yr“\ which we claim as the true mass loss rate for 77 Carinae for most of the 5.52 year cycle. 
However, the mass loss rate from 77 Carinae is not well constrained by observations in other regions 
of the EM spectrum. Andriesse, Donn and Viotti (1978) derived a value of M,, ss 7.5 X 1O~ 2 M 0 yr -1 
based on their analysis of the dust formation rate, while van Genderen and The (1984) estimated 
that M n ss 1 O - 4 AT 0 yr -1 based on a number of arguments. Davidson (1987) gave M n < 1O _ 3 M 0 
yr -1 . Radio observations (White et al. 1994) yield M v » 3 x 1 O _ 4 M 0 yr -1 . A STIS spectrum 
obtained in March 1998 (coincidently near the time of our “enhanced mass loss” interval) and 
analyzed by Hillier et al. (2000) yielded a value of M n ss lO“ 3 A / 0 yr" 1 . Millimeter observations 
(Cox et al. 1995; Cox 1997) yielded values of M ^ ss 1O _ 3 M 0 yr -1 . We should point out that most 
of these other rfiass loss rate estimates assume that 77 Carinae is a single star; this assumption 
neglects the effect of the companion on the ionization balance in the outer wind of 77 Carinae which 
in turn can significantly affect mass loss estimates which may be based on the strength of emission 
line profiles or the amount of free-free emission in the outer part of the wind. It, is possible that 
by including the influence of the companion on the ionization of 77 Carinae’s wind, and/or the 
inclusion of wind variability, better agreement in the values of mass loss rate derived from analyses 
in different wavebands may be obtained. 
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While wo have modelled the X-ray eclipse as an interval of enhanced X-ray absorption, an 
alternative possibility is that the long X-ray minimum might represent a reduction of the intrinsic 
X-ray emission. Such a reduction presumably could be produced if the winds near periastron collide 
at velocities much less than their terminal speeds. The relative collision velocities could be less 
near periastron if the winds from one or both stars have slower accelerations than winds from other 
massive stars, or if radiative braking (Gayley, Owocki, and Cranmer 1997) is important. 

We expect that STIS (or other high spatial resolution) Paschen-line observations of the core 
of 7 / Carinae near the next periastron passage (2003.5) will show large velocity shifts if the binary 
model presented here is correct, and the mass ratio is not too extreme. However, if significant 
variations in the ionization, the density and/or emitting volume near periastron passage occur, 
then determining accurate radial velocities and orbital elements from the emission lines may be 
difficult even with STIS. The X-ray fluxes however should be relatively insensitive to these effects. 
It may very well be that the variable P-Cygni absorption components seen in the STIS data (Hillier 
et al. 2000) may provide a good estimate of the stellar velocities, since the absorbing volume is 
constrained to lie in front of the stellar disk. Coordinated X-ray and STIS spectral observations of 
the P Cygni lines near periastron will be especially useful to fully constrain the orbit. 
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Fig. 1. Plot of RATE count rates versus contemporaneous ASCA 2-10 keV source flux. The 
RXTE rates have been corrected for instrumental background. The straight line shows the best fit 
to the data. 


Fig. 2.— The RXTE lightcurve of ?? Carinae and the colliding wind model (thick solid line). The 
thick dotted line is the model derived using the orbital elements in Damineli et al. (2000), with 
f(m) — 7.5. The model based on the Damineli et al. (2000) parameters does not fit the abrupt 
decline in the observed 2-10 keV X-ray flux. Neither model fits the width of the X-ray eclipse 
adequately. The thin dotted line shows the effect of increasing M n to 3 X 1 0~' f Mq yr _1 . 

Fig. 3. The RXTE lightcurve of i) Carinae and the colliding wind model (see Table 1), including 
an interval of enhanced mass loss following periastron passage using the parameters in Table 1, 
along with an enhanced mass loss model based on the Damineli et al. (2000) parameters. The extra 
mass loss from the primary has the effect of suppressing the observed flux from the system for an 
extended period following periastron passage, and provides a better match to the duration of the 
X-ray eclipse. 

Fig. 4. The variation in TV// to the colliding wind shock. The “enhanced-M^” model is shown by 
the thick line, and the “constant- M v " model is shown by the thin line. The “enhanced- A/,/’ model 
provides an extended interval of extra absorption (N H > 10 24 cm -2 ) which completely hides the 
colliding wind shock. 

Fig. 5. The expected radial velocity curve of the primary based on the orbital elements derived 
in the present work (solid line) and on the Damineli et al. (2000) orbital elements (dotted lino). 
The ?> 0 .7 velocities from the STIS data reported by Davidson et al. (2000) are also shown. 



RXTE Net Count Rate (3 PCUs) 






12 






Radial Velocity (km s’ ) N H (cm 











- 15 - 


Table 1. Colliding Wind Model Parameters 


Parameter 

Present Work 

Damincli et al. 2000 Value 

T (periastron) 

1997.95±0.05 

1998.13 

e 

0.90 

0.75 

P 

5.52 years 

5.53±0.01 

Mr, 

80 M© 

70 

¥c 

30 Mq 

68 

to 


275° 
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-12 km s 1 

Mr, 


i— 
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Me 


10- 5 M© yr- 1 

V 


500 km s _l 

v 

v oo ,c 


2000.0 km s’ 1 



